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We study mean value properties of harmonic functions in metric measure spaces. 
The metric measure spaces we consider have a doubhng measure and support a (1,1)- 
Poincare inequahty. The notion of harmonicity is based on the Dirichlet form defined 
C^ ' in terms of a Cheeger differentiable structure. By studying fine properties of the 

Green function on balls, we characterize harmonic functions in terms of a mean value 
property. As a consequence, we obtain a detailed description of Poisson kernels. We 
shall also obtain a Gauss-Green type formula for sets of finite perimeter which posses 
a Minkowski content characterization of the perimeter. For the Gauss-Green formula 
we introduce a suitable notion of the interior normal trace of a regular ball. 
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1 Introduction 

Solving the Dirichlet problem on a smooth domain in R" is equivalent to constructing har- 
monic measure on the boundary of the domain. More precisely, it is known that the classical 
harmonic measure can be expressed in terms of a Poisson kernel which is given by the Radon- 
Nikodym derivative of harmonic measure with respect to the Hausdorff boundary measure; 
that is 

In general metric measure spaces with a doubling measure and a Poincare inequality the 
Dirichlet problem has been solved for Sobolev type boundary data in |Sh2] . and also for all 
continuous boundary values in |BBSj . In |BBS] the authors provide an integral representation 
for the solution to the Dirichlet problem, and hence extend the solvability to L^ boundary 
data. In this general setting, however, due to lack of a natural choice of boundary Hausdorff 
measure one has to replace a Poisson kernel with a Poisson kernel-like object for which 

Pxo{x,y) = -7—{y)- 

It was shown in [BBSj that for a fixed Xq G Q, where fi is a bounded open subset of X, 
there exists a Radon measure Uxg concentrated on dQ, i.e. Uxg is a harmonic measure on dQ 
evaluated at xq, and a real- valued function P^^ on f2 x dfl such that whenever / G L^{dQ, u^q) 
the following expression for the harmonic extension Hj is valid: 



Hfix) = f(y)Pxo{x,y)dux,Xy), 
Jan 



Ida, 

and moreover, for each y G dQ the function P^gi'^y) is harmonic in Q. 

Our main objective is to find a relationship between the Poisson kernel that generates 
solutions to the Dirichlet problem in terms of Cheeger differentiable structure, and the 
perimeter measure of a ball of finite perimeter in metric measure spaces. Our framework is 
a complete geodesic metric measure space with a doubling Borel measure, and we moreover 
assume that the space supports a (1, l)-Poincare inequality. These conditions are discussed 
in detail in Section [2l We shall describe the Poisson kernel in terms of an analog of a normal 
derivative of the Green function at the boundary. 

We also study divergence-measure fields along the lines of Ziemer [Z] in this general 
context. We consider an L^-vector field, F, from a metric measure space X to M'^ for which 
div F is a real- valued signed Borel measure with finite mass. 

To investigate divergence-measure fields we shall provide a meaningful definition for the 
divergence operator in metric measure spaces. We then generalize some results obtained in 
[Z] to the metric setting. In particular, we obtain the Gauss-Green type integration by parts 
formula for sets of finite perimeter which possess a Minkowski content characterization of 
the perimeter. For the Gauss-Green formula we introduce a suitable notion of the interior 
normal trace of a regular ball. 

We mention a related paper by Thompson and Thompson |TT] in which the authors 
define divergence and prove an analogue of the Gauss-Green theorem in Minkowski spaces, 
i.e. in finite-dimensional real normed spaces with smooth and strictly convex unit ball. 



We use the results for the divergence operator to characterize the Laplace operator of 
the Green function on regular balls as the sum of the Dirac point mass and a measure 
concentrated on the boundary of the ball. This characterization allows us to give a precise 
description of the Poisson kernel defined in |BBS] . In the setting of Heisenberg groups, we 
explain the relation between this measure and the perimeter measure or the codimension 
one Hausdorff measure. 

2 Preliminaries 

Here we recall some basic definitions and the notation we shall use in this paper. Our 
framework is given by a complete metric measure space {X,d,fi), where fi is doubling, that 
is, there is a constant c > such that for every ball B = Br{x), x E X and r > 0, 

< fx{2B) < cfi{B) < cx). (1) 

We write Br{x) for the ball centered at x with radius r > 0, and XB = Bxr{x) for any A > 0. 
The smallest value of c for which ([1]) is valid is called the doubling constant of X, and we 
shall denote it as q. 

An upper gradient for an extended real-valued function m : X — )■ [— oo, +oo] is a Borel 
function g : X ^ [0, C)o] such that 

|M(7(0))-M(7(g)|< [gds (2) 

for every nonconstant compact rectifiable curve 7 : [0, l^] — )■ X. We say that g is a. p-weak 
upper gradient of u if ([2]) holds for p-almost every curve; the notion of p-almost every curve 
is in the sense of the p-modulus of a curve family F defined as 

Modp(r) = inf < / g'' dfi : g > is a Borel function, / f? ds > 1 for all 7 G F > . 

If u has an upper gradient in U'{X, /i), then it is possible to prove the existence of a unique 
minimal p-weak upper gradient gu G LP{X,fi) of u, where gu "^ g iJ'-a.e. for every p-weak 
upper gradient g G L'p{X,jj) of u. We refer to |Sh2] for the case p > 1, and for the case 
p = 1 to |Haj03| . 

In what follows, the metric space is supposed to support a weak (1, l)-Poincare inequality: 
there exist constants c > and A > 1 such that for all balls B^- with B\r C X, for any 
Lipschitz function / G Lip(X) and minimal p-weak upper gradient gj of / we have 

J \f ~ fBr\ dfx<cr j- gf d/x, (3) 

J Br J B\r 

where 



1 

'Br fJ-iBr) JB, 

is the integral average of / on Br{x). 



fBr ■= 4 fd^i := -j^ I fdfi 

J B, 



It is well known that the doubling condition and the Poincare inequality imply the quasi- 
convexity of the metric space X, see [K] and |HaKo] . Therefore, up to a bi-Lipschitz change 
of the metric, the space X can be assumed to be geodesic, that is, given x,y E X there is 
a curve 7 with end points x, y and length d{x,y). Moreover, for a geodesic space the weak 
(1, l)-Poincare inequality implies the (1, l)-Poincare inequality, i.e. ([3]) holds with A = 1. 
Therefore, as most of the properties of metric spaces we consider are bi-Lipschitz invariant, 
it is not restrictive to assume that X is a geodesic space and supports a {1,1) -Poincare 
inequality. 

As proved by Cheeger in [C] , in our setting the following differentiable structure is given. 
There exists a countable measurable covering Ua of X, and Lipschitz coordinate charts 
X" = (Xf , . . . ,X^J : X ^ M'^- such that //([/<,) > for each a, /x(X \ [j^ U^) = and 
for all a the following holds: the charts (X°, . . . ,X^ ) are linearly independent on 11^ and 
1 < ^a < X, where X is a constant depending on the doubling constant and the constants 
from the (1, l)-Poincare inequality satisfying the following condition: For any Lipschitz 
function / : X — )■ M there is an associated unique (up to a set of zero /i-measure) measurable 
function daf '■ Ua -^ ^^°' for which the following Taylor-type approximation 

f{x) = /(xo) + daf{xo) ■ (Xaix) - Xc,(xo)) + o(d(x, Xo)) (4) 

holds for /i-a.e. Xq E Ua- 

The previous construction implies, in particular, that for x E Ua there exists a norm \\-\\x 
on M*^" equivalent to the Euclidean norm | • |, such that gf{x) = \\daf{x)\\x for almost every 
X E Ua- Moreover, it is possible to show that there exists a constant c > 1 such that 

-gf{x) <\df{x)\< cgf{x) 

for all Lipschitz functions / and fi-a-e. x G X. By df{x) we mean daf{x) whenever x G Ua- 
Indeed, one can choose the cover such that Ua H Uj3 is empty whenever a ^ (3- 

Formula (jlj) implies in particular linearity of the operator f \-^ df and also the Leibniz 
rule d{fg) = fdg + gdf holds for all Lipschitz functions / and g. 

For the definition of the Sobolev spaces N^'^{X,^) we will follow [ShlJ. Since we assume 
X to satisfy the (1, l)-Poincare inequality, the Sobolev space X^'^(X, /x), 1 < p < 00, can 
also be defined as the closure of the collection of Lipschitz functions on X in the following 

jVi.p_norm 

II IIP II IIP I II IIP 

\\^\\i,p — \\^\\lv{x) ^ lly«llLP(x)- 

The space X^'^(X, /i) equipped with the X^'^-norm is a Banach space and a lattice [Shlj . 
Let ii^ C X be a Borel set. The p-capacity of E is defined as usual to be the number 



Capp(E)=inf ( / |M|Prf/i+ / \du\Pdii 

where the infimum is taken over all u G X^'P(X, yu) for which u = 1 on E. We say that a 
property holds p-quasieverywhere, p-q.e. for short, if the set of points for which the property 
does not hold has p-capacity zero. For instance, if 'U,f G N^'^{X,fi) and u = v /x-a.e., then 
u = V p-q.e. and ||-u — f ||i,p = 0. If we, moreover, redefine a function u G X^'P(X, /x) on 



a set of zero p-capacity, then it remains a representative of the same equivalence class in 

We shall also use Sobolev spaces defined on a domain Q (i.e. a non-empty open path- 
connected set) of X; the space N^''^{Q, /i) is defined in the same way the space A^^'^(X, /i) is, 
but considering Q as the ambient space. The space of Sobolev functions with zero boundary 
values is instead defined as 

N^'^'in, fi) = {u G N^'P{X, /i) : u = p-q.e. on X \ (]} . 

We have that Nq'^{Q, fi) = N^''p{X, ^) as Banach spaces if and only if Capp(X \ fi) = 0. 

In what follows, let p = 2. By [FHKj . the Cheeger differentiable structure extends to all 
functions in N^''^{X,fi) and N^''^{Q, fi), and hence we define an inner product on N^''^{X,fi) 
by the Dirichlet form 

£{u,v) = I {du,dv)dfi, 
Jx 

for all M,f G N^''^{X, fi). It can be proved that such a form is strongly regular with the 
domain, or core, given by N^''^{X, fi). 

We recall that a Dirichlet form is said to be strongly regular if there exists a subset K 
of the domain of the Dirichlet form, dense in both this domain and in the class of Lipschitz 
functions on X, such that the distance dg : X x X ^ [0, oo] defined, in our case, by 

deix, y) = sup {ip{x) - ip{y) : \dip{x) | < 1} 

is a metric on X that induces the same topology on X as the original metric topology on X. 
In fact, under the doubling property and a Poincare inequahty dg is bi-Lipschitz equivalent 
to the original metric d on X, and so the Dirichlet form 8{u, v) is strongly regular. The set 
K is called a core of £. We refer to |Sturm] and |FOT] for more details. 
For each a > we define the bilinear form 



£a{u,v) = a uvdn + £{u,v). 
Jx 

We thus have on N^'^{X,ii) the norm || ■ \\a induced by Sa which is equivalent to the N^'^- 
norm. In this way, N^'^{X,fi) with the norm || ■ \\a is a Hilbert space with inner product Sa- 
Note that S by itself is not an inner product on N^''^{X, /i); S{u, m) = if and only if m is a 
constant (see [C]). If, for example, fJ,{X) < oo, then £{u,u) = does not imply that u = 0. 
The fact that the bilinear form S^ yields a Hilbert space can be seen as follows. Since 
the X^'^-norm is comparable to the £^Q-norm, we have that N^''^{X,fi) is complete also with 
respect to the £^Q-norm. In this way the ^c^-norm is well defined for any u G X^'^(X, /i). By 
approximation and the linearity of the map u i— >■ du, the Leibniz rule follows for functions u 
and V in N^''^{X,fi) (we refer for these properties to the paper |FHKj ). 

Remark 2.1. We point out that the convergence of a sequence {uk)k to a function u in 
N^''^{X,^) is same as the convergence of the two sequences {uk — u)k and {guk-u)k to in 
L2(X,/i). 

In general, the convergence of Uk to u in L'^{X,fi) together with the convergence of gu,. 
to gu in L'^{X, /i) does not imply that Uk converges to u in N^''^{X, /i). As a counterexample. 



consider the metric space X = M^ with the distance induced by the norm ||(x,y)||i = |x| + ||/| 
and with fi the Lebesgue measure; in this case the upper gradient is determined by the dual 
norm ||(a;,?/)||oo = max{|a;|, \y\}. It suffices to verify this for a Lipschitz function u. For such 
function, by [C], denoting by Br {xo,yo) the ball in the norm || ■ ||i with radius r centered 
at (0:0,^0)7 "we have that 

. ^ ,. \uix,y)-u{xo,yo)\ 
gu{xo,yo) = hra sup 

{x,y)&B). '(xo,yo) 

= max (Vm(xo,i/o),^^)r2 = ||VM(a;o,i/o)||oo- 

||^||i=i 

The sequence Uk{x,y) = x + fk{y), where fk{y) = dist(?/, ^Z) converges to the function 
u{x,y) = X, but for a.e. point 

9uAx,y) = \\Vuk{x,y)\\oo = 1= \\\/u{x,y)\\^ = gu{x,y) 

and 

gu^-uix,y) = \\\/ukix,y) -Vu{x,y)\\oo = ||V/fc(y)||oo = 1- 

Nevertheless, it is possible to use Mazur's lemma to prove that for a convex combination 
the aforementioned property holds true, both for the Cheeger differentiable structure and 
for the upper gradient. For the Cheeger differentiable structure, however, it is not necessary 
to take convex combinations. Indeed, in this case the sequence of gradients duh is bounded 
in L^(X, M'^,/!), and so it is weakly convergent to some if G L'^{X,W^, fi). Mazur's lemma is 
then needed only to show that ip = du. We can consider convex combinations 

N{k) 
i=l 

with strong convergence Vk ^ urn. L'^(X, /i) and dvk -^ (p in L'^{X, M.'', ft), that is Vk ^ u in 
N^''^{X, fi), and we may then conclude that ip = du. We then obtain 

lim / \du — duk\'^ djj, = I \du\'^diJi+ lim I / \duk\'^ dfj, — 2 I {du,duk)dn 
^^°° Jx Jx ''^'^ \Jx Jx 



2 / \du\^diJ,-2 / {du,ip)dij = 
'x Jx 



by the weak convergence. 



3 Metric Laplace operator 

In this section we construct a metric Laplace operator A^ on the metric measure space 
{X, d, fi). Recall that a Dirichlet form S is strongly local if whenever u, v are in the domain 
of S and u is constant on the support of v, then S{u, v) = 0. Having a strongly local Dirichlet 
form at one's disposal it is rather standard argument to construct an operator associated 
to the form. Most of the statements (without detailed proofs) can be found in the book 

6 



of Fukushima, Oshima and Takeda |FOTj . but we provide complete proofs for the reader's 
convenience. Since this operator plays the role of the Laplace operator on X, we shall denote 
it by Ax- Setting 

Dom(Ax) = {ue N^''^{X,fi) : there exists / G L'^{X,fi) 

such that £{u,v) = — / fvdfi for all v G A^^'^(X, yu)} , 

Jx 

the Laplace operator is defined by 

Axu = f. 

We summarize the main properties of this operator in the following theorem. The main 
point is to construct the resolvent operator Ra, i.e. an operator that gives for any a > the 
formal solution of the problem 

(a - Ax)u = /, (5) 

and to deduce from this the main properties of A^. 

Theorem 3.1. For each a > 0, there is an injective bounded linear operator Ra : L'^{X, fi) — )■ 
iV^'^(X, n) such that for all v G N'^''^{X, /x) 

/ fvdn = Sa{Raf, V) = ^{Raf, V) + a{Raf, v)2- 

Jx 
This operator satisfies: 

1. foranyfEL^iX,fi), WRJh < ^\\fh; 

2. for any a,P > 0, Ra{L^{X,fi)) = Rp{LP'{X^ ^)), and the resolvent equation holds true 

Raf - Rpf = (/3 - a)RaRpf (6) 

for all f G L^{X,fi); 

3. for any f G L^(X, yu), we have the following limit in the L'^{X, fi)-norm; 

lim aRaf = f (7) 

Properties 2. and 3. imply that Ra{L^{X,^)) is dense in L^(X, /i). In addition, 

Dom{Ax) = Ra{L\X,fi)) 
for any a > 0, and for u = Raf , AxU := au — f is independent of a. 

Proof. Let us fix / G L^{X,fi). Then we can define the linear operator Tf : A^^'^(X) — )■ M 
by Tf{v) = (/, v)2 := Jj^ f V dfi. We have that 



\Tfiv)\<\\fh\\vh<^Sa{v,Vy/'. 

\/Oi 



Therefore Tf is a bounded linear operator on the Hilbert space {N^''^{X, ^),£a), so by the 
Riesz representation theorem, there exists an element of A^^'^(X, /i), denoted by Raf, such 
that Tf{v) = £a{Raf,v). The map Ra '■ L'^{X,^) — )■ A^^'^(X, /i) defined above is linear by 
the linearity of the defining operator f \-^ Tf. 
Since 

a{Raf, V)2 = £a{Raf, v) - £{Raf, v) = (/, v)2 - £{Raf, v) , 

choosing v = Raf and applying Holder's inequality, we see that 

< a\\RJ\\l = a{RJ, RJh = if, RJ)2 - S{RJ, RJ) < (/, RJ)2 < WfURJU 
Thus we obtain Claim 1 of the theorem, namely, 

a\\Rafh<\\f\\2- 

Thus Ra as an operator mapping L^{X, /i) to L'^{X, fj,) is bounded with image in N^''^{X, fi) C 
L'^{X,fi) and its operator norm given by 

\\Ra\\ '■= WRoWl'^^L^ < — ■ (8) 

a 

We now prove the resolvent equation ([6]). Let us take / G L^{X,fj,) and v G N^''^{X, fi). 
Then 

£a{Raf - Rpf +{a- /3)RaRpf, V) = £a{Raf, v) - £a{Rf}f, v) + (a - /3)£a{RaR^f, v) 

=(/, vh - £{Rpf, v) - a{Rpf, vh + {a- P){Rpf, v)^ 
=if,v)2-£p{Rfif,v) = 0. 

This means that for / G N^''^{X, fi), and then by density also for / G L'^{X,fi), we have the 
identity 

Raf - Rpf + (a - m^Rpf = 0. 
Moreover, if we consider / G A^^'^(X, /x), we have (denoting £a{f-, fY^"^ 

a\\aRaf - f\\l< £a{aRaf - /, aRaf - f) 

= a^£a{Raf, Raf) + \\f\\l ' '2a£a{Raf, f) 

= a\f,Raf)2 + £{f,f)-a\ 



By Holder's inequahty and by using ([8]), we also get 

a{f,Raf)2 - WfWl < C^WRafh ll/lh - 11/11^ < 0. (9) 

Therefore, 

lim WaRaf - fh < lim J^iMl = q. 

a— >oo a— >oo y tt 

To extend this limit to be valid for any / G L^(X, /x), we use the boundedness of Ra, by 
fixing /j G A^^'^(X, /i) such that ||/ — /£||2 < £^- In this way we get that 

WaRaf - fh < Wc^Rafe ' fsh + «Pa(/ - fe)h + 11/ " feh < \\aRafe ' feh + 2^, 



and hence 

limsup \\aRaf — /II2 < Se. 

From this Claim 3 of the theorem follows since e was arbitrary. 

We have now proved that Ra is a strongly continuous resolvant (see [FOT] ) for any a > 0. 
Let us next prove injectivity of Ra- Suppose / G L^(X, /x) is such that Rjsf = for some 
(3 > 0. Then by the resolvant equation ([6]), 

= Raf- Rpf + (a - l3)RaRpf = RJ. 

that is, Roif = for every a > 0. Now by equation ([7]), we see that / = 0, that is, Ra 
is injective. We can therefore define the inverse map R~^ : Rai^L"^ {X , fi)) — )■ L^(X, /i). We 
claim that 

Dom(Ax) = Ra{L^{X, /i)), A^m = au- Ra^u. 

For this definition to be consistent, we first show that the set Ra{L'^{X, /i)) and the operator 
AaU := au — R'^u do not depend on a. By the resolvent equation ([6]), 

Rpf = RaU+{o^-mpf)- 

Therefore, for every / G L'^{X,fi), R/sf G i?Q,(L^(X, /i)), and hence 

Rp{L\X,fi))cRa{L\X,fi)). 

By the symmetry of the argument, we have the required result Ra{L'^{X, /i)) = Ri^^L'^^X, yu)). 
Let us write D = Rp{L'^{X, ^)). 

li u E D and a, (3 > 0, then AaU — Apu = {a — (3)u — R~^u + R'^^u. Therefore, 

Ra{AaU — Ai3u) = aRaU — (3RaU — U + RaR^^U. 

On the other hand, since D = i?^(L^(X, /i)), there exists / G L^(X, /i) such that i?^/ = u. 
Hence we have, by the resolvent equation ([6]), that 

Ra{AaU - Apu) = aRaRpf - PRaRpf - R/sf + Raf 
= Raf - R(sf + (« - P)RaRpf = 0. 

By injectivity of Ra, we see that AaU — Apu = 0, i.e. AaU = Apu. 

Let us now show that i?Q,(L^(X, /i)) C Dom(Ax). Let u G Ra{L'^{X, ft)). Then there 
exists / G L'^{X,fi) such that u = Raf- The identity £a{Raf,v) = (/, f)2 for any v G 
N^''^{X,fi) can be written out as follows 

{du, dv) dfi =S{u, v) = £a{u, v) — a(u, v)2 

X 

= £a{Raf,v) -a{u,v)2 

= (/, v)2 - a{u, v)2 = - / {au - f)v djj 

Jx 



for all V G N^''^{X, ^). This simply means that u G Dom(Ax) and that Axu = au — f - 
au — R~^u = AaU. 

For the reverse inclusion, Dom(Ax) C Ra{L'^{X, fi)), let us consider u G Dom(Ax; 
Thus there exists / G L'^{X,fi) such that for all v G N^''^{X,fi) we have 

{du, dv) dfi = — fv dfi. 
X Jx 

Then consider w := Ra{au — /); we obtain that 



£a{w,v) = {au — f,v)2 = a uvdfi— / fvdfi = a uv dfi + £{u,v) = £a{u,v), 

Jx Jx Jx 

that is w = u, which means that u G -Rq,(L^(X, /x)). The identity / = Axu = au — R~^u 
follows easily. 

In addition to the density of Dom(Ax) in L'^{X, ^), we also have that Dom(Ax) is dense 
in A^i'2(X,/i). In fact, by ® and ([9]), for any / G N^'^{X,ij,), we have that 

Wc^Raf - fWl =£a{aRJ - /, aRJ - /) = a^£^{RJ, RJ) - 2a£^{RJ, f) + £M, f) 

= «'(/, Raf)2 - «(/, /)2 + £{f, f) < £{f, /), 

that is the sequence {aRaf — f)a is bounded in N^'^{X, fi). Therefore, for any sequence of 
positive real numbers (a„)„ so that hm„a„ = cxd, the corresponding sequence of functions 
C(nRa„f — / is a bounded sequence, and hence by Mazur's lemma we have a sequence of 
convex combinations converging in N^''^{X, fi); 

'N{n) \ 

J=n ) 

On the other hand, limQ_^oo «-Ra/ = / in L'^{X,f^)- Thus w = /i-a.e. in X, and hence 
by the fact that w G N^'^{X,fi) we know that w = p-q.e. in X. Therefore, it must be 
that w = 0. Observe that the sequence of convex combinations ^jj„ Xi^nO^iRaJ lies in 
Dom(Ax) and converges to / G A^^'^(X, /i), so the proof is completed. D 

We can now give the definition of a Cheeger harmonic function in the obvious way. 

Definition 3.2. A function u G A^^'^(f2,/x) is said to be Cheeger harmonic (referred to in 
this paper as harmonic) if 

{du, dv) dfi = 



rl,2/ 



for all V G Nq' (fi, /i), i.e. u is harmonic if and only if m G Dom(AQ) and Aqu = 0. Here A^ 
is the operator defined in Remark 13.41 below. 

Remark 3.3. The notion of Cheeger harmonicity refers to the fact that we are using the 
Cheeger differentiable structure. This notion has been previously considered in the paper 
[KRS] , where Lipschitz regularity of Cheeger harmonic functions has been investigated. We 

10 



also underline that Cheeger liarmonicity can be equivalently be given in terms of a minimizer 
of the Dirichlet energy: u is Cheeger harmonic if and only if for any ball B,. 

\du\^ d^ < I \dv\^ d/i, 

for all V such that v — u ^ Nq' {Br, /i). 

Remark 3.4. Let Q G X he a. bounded domain satisfying a (1, 2)-Poincare inequality with 
Cap2(X \ f]) > 0. The previous construction of Ax can also be used to construct a Laplace 
operator on the subdomain Q. There are essentially two different Laplace operators; the first 
is just the restriction of Ax to Q and is defined by 

Dom(An) = {u E N'^''^{Q,,^) : there exists / G L^(f2,/i) such that 

/ {du, dv) d^= — I fvdfi for all v E Nq' (fi, /i) > , 
Jn Jn ) 

and the operator is given by 

Ann = /. 

The second alternative, adapted to the inhomogeneous Dirichlet problem, is the operator 
defined by 

Dom(A^) = {m G A^o' (^)/^) • there exists / G L^(f2,/i) such that 

/ {du, dv) dfi = — / fvdfi for all v E N^' (fi, fi) > , 
Jn Jn } 

and 

Ku = f. 

To define the latter operator, the previous procedure has to be modified by considering the 
Hilbert space A^q' {fi,fj.) with the inner product £a for all a > 0, to obtain the resolvent 
operator i?° : NQ'^{n,fj.) -^ iVo'^(n,/i) with 

S^{Rlf,v) = {f,vh 

whenever v G Nq' {fl,fi). Since the vector subspace Lipo(fi) of Nq' (fi,/i) is also a dense 
subspace of L'^{fl,fi), we may extend i?° to be an injective map from L'^{Q) to Nq' {Q). 
These properties, like the one proved for R^ in Theorem I3.H are the essential properties for 
the definition of the operator A^. 

It is easy to verify that the operator Aj^ is the restriction of Ax to Q in the following 
sense: If m G Dom(Ax), then 

^n{u\n) = (Axm)|q. 

On the other hand, the operator A^ is the restriction of A^ to the space A^g' {Q, /i), that is 

Dom(A^) = Dom{An)nN'Q'\Q,fi) 
with A^u = Anu for u G Dom(An) n iVo'^(fi,/i). 
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3.1 Measure- valued Laplace operator 

Let fi be a domain in X. We give the following definition of the measure- valued Laplace 
operator ^^ on Q. By ^i){Q) we denote the space of all bounded signed Borel measures on 
Q, i.e. z/ G ^b{^) is a real-valued signed Borel measure on Q with bounded total variation 



|z/|(n) = sup < / (pdu : (p E Lip^ 



(fi), ||^||oo<U <00. 



We remark that to compute the total variation of a measure we test in the space Lip^(r2) 
of Lipschitz functions on Q with compact support instead of the space Cc{^) of continuous 
functions with compact support; we may do this since Lip^(r2) is clearly dense in Cc{fi). 
We define 

Dom(^n) ={u e N^''^{Q,fi) : there exists u G ^h{X) such that 

£(u, v) = - I vdv for all v G Lip^(fi) I , (10) 

and then we set 

^qu = z/. 

Example 3.5. As an example, we can consider the Euclidean space {MJ^, || ■ ||) and modify its 
metric structure in two ways, which essentially lead to the same metric measure structure. 
We fix r2 C M" an open set with regular boundary and we can modify either the measure by 
considering dfi = {1 + xn)dC"', or the differential structure du = {1 + xn)VM. In both cases 
we have for u,v G C^(]R") 

{du, dv) dfi= Vm ■ Vf dx + Vn ■ Vf dx 

Jr" Jn 

vAudx— / vAudx+ / vVu ■ undT-L"^^. 
Jn Jan 

Then u G Dom(A]Rn) if and only if Vu ■ z/f^ = on dft and Am G L^(R"'). In addition, in the 
case /i = (1 + Xn)^^ with the standard differential structure we also have 

AiRnM = Am. 

In the second case, where fj, = C"^ and du = {1 + Xn)'^u, we obtain 

A^nu = (1 + Xn)Au. 

In a similar fashion, Dom(^]Rn) is given by those functions u for which Au G L^(]R") and 
the trace of Vu ■ vq G 1.^(9^2,7/""^), and 

^MnM = Aw/i - Vm ■ z/n?^"-^L9fi. 

It can be verified that Dom(^o) is a vector space and that ^n is linear. We wish to 
expand the class of test functions in the definition of the domain Dom(^Q) from Lip^(r2) to 
allow for test-functions t> in A^q' (f2,/x), see Proposition 13.81 For that, we need the following 
lemma. 
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Lemma 3.6. If E G Q is a Borel set such that Cap2(-E) = 0, then for every u G Dom(^Q), 

Proof. By the Jordan decomposition theorem, the measure S>nu can be decomposed into its 
positive and negative parts, S^^u and S>^u\ this means that we can decompose Vt into two 
disjoint Borel sets Vt = Q^ U Q^ in such a way that ^qu{B) > for every B C Q^ and 
^qu{B) < for every B C O^. Hence we may, without loss of generality, consider E C fi^; 
in fact we can decompose E = E^ U E~ and use the monotonicity of capacity. Further, we 
may also assume that E' is a compact set, since as Radon measures both 2>^u and S>^u are 
inner measures and E is a. Borel set. 

Since Cap2(-E') = 0, we have also that the relative capacity Cap2(-E',fi) is zero. This can 
be seen by multiplying those Lipschitz test-functions which were used for computing Cap2(-E') 
by another Lipschitz function rj which is 1 on a neighborhood of the compact set E and has 
compact support in Vt. We can then find a sequence of Lipschitz functions {ipi)i so that 
< (/'i < 1 on X, (/9j = 1 on i?, and ||(y9j||jvi.2(x) < 2~*, and ipi are compactly supported in Vt. 
We may assume that the sequence {^i)i converges pointwise to zero outside of the compact 
set E (we can do so by choosing ipi to have support in the open set IJa;e-B-^('''' lA))- We 
have 



X 



1/2 / ^ \ 1/2 

< ( / Idul'^djj) / \d(pi\'^dn 



X J \JX 



{du, difi) dfi 

'X 
< I / \du\^dfi\ ||v5i||Ari.2(x,M), 

which tends to as z — ?■ oo. 

On the other hand, since (pi are all bounded by 1 and |^n|(-^) < oo, by the Lebesgue 
dominated convergence theorem we have 



lim / ipid^nu = ^nu{E) = &^u{E). 



A similar argument shows that ^^u{E) = 0, and hence the proof follows. D 

Remark 3.7. The requirement that E is a Borel set in the above lemma is not a serious 
restriction, because ii E G Q is a set with Cap2(-E') = 0, then there is a Borel set £"0 with 
E C Eq cQ such that Cap2(-Eo) = 0. 

The following proposition tells us that we do not have to restrict ourselves to having 
test-functions v only in Lip^(X) in f lTO|) . 

Proposition 3.8. Letu G Dom(i^f7). Then for every v G A^q' {Q, fi)r\L°°{Q, ^) the following 
holds: 

£{u,v) = — V d^QU. 
Jq 

Proof. We first assume that v ha compact support in Q. Note that by the (1, 2)-Poincare 
inequahty we can approximate compactly supported functions in A^^'^(f2,/i) by Lipschitz 
functions. So we can find a sequence of compactly supported Lipschitz functions {ipi)i on Q 
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that converge to v in the N^''^{VL, /i)-norni. By passing to a subsequence if necessary, we may 
also assume that ipi ^ v pointwise outside a set of zero 2-capacity; we refer to [Shlj . Since 
V is bounded, we can also assume that the approximating compactly supported Lipschitz 
functions ipi are also uniformly bounded by M := ||f ||oo- Applying ipi as in (ITOl) . we see that 



(Pid^QU = — / {dipi, du) dfi ^ — / {dv,du) d^. 
n Jq Jn 

By Lemma 13.61 we know that ipi —^ v almost everywhere with respect to the total variation 
measure |^qm|. By the Lebesgue dominated convergence theorem applied to the uniformly 
bounded functions yjj with respect to the positive and negative parts ^qm"*", 2i^u~ of the 
signed Borel measure S)^u^ we may conclude that 



n Jn 

Hence equation fITOl) holds for all compactly supported functions v G N^'^{Q,fi) H L°°{Q, fi). 
To pass to any v E N^' (f2,/i) fl L°°{Q,fi), we note that functions in A^q' (f2,/x) with 
compact support in Q form a dense subclass of Nq' (fi,/i) (see |Sh2] ). Hence, if v is in 
Nq (fi,/u) n L°°(fi,/i), we can find a sequence of compactly supported functions Vi from 
A^q' (fi,/u) n L°°{Q,fi) such that Vi -^ v in Nq' (fi,/i). As before, we can also ensure that 
Vi —^ V 2-capacity almost everywhere in Q. Hence 

{du, dv) d^ = lim / {du, dvi) d^ = — lim / Vi d!^nu, 
n *-^°° Jn *^°° Jq 

and then if v is bounded in Q we have 

lim / Vi d!^QU = v d!^QU, 
»-^°° Jn Jn 

giving the desired result for all bounded functions in A^g' (fi,/i). D 

We shall also need the following lemma, which is based on the Lebesgue decomposition 
of the measure ^n given by 

d^nu = fudfx + d^^u, 

where /„ = ^^^^ is the absolutely continuous part and ^q the singular part of !^q. 

Lemma 3.9. Let u G Dom(^n). // the singular part S^^ of S>q,u is zero and if the Radon- 
Nikodym derivative /„ G L^{Q,fi), then u G Dom(An) with Aqu = fu- 

Proof. From the discussion in Section |3l if ^^u = and the absolutely continuous part is 
represented by /„ G L^{Q,fi), then 

S{u,v) = - / fuvdfi (11) 

Jn 

for all V G Nq' {Q, fi)nL'^{Q, n). When /„ G L'^{Q, fj,), we can use a truncation argument and 
the Lebesgue dominated convergence theorem to show that (ITT]) holds for any v E Nq' [fl, fj,). 
So we conclude that u G Dom(Ao) and A^u = fu- 

n 
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Remark 3.10. It can be seen that Doni(Af7) C Doni(^n); moreover, ii u,v G Dom(^Q) 
and a G M, the following hold true: 

1. spt(^Qu) C spt(u); also, if u is constant on an open set U, then spt{^nu) C Q\U; 

2. u + v,aue Dom(^n) with ^n{u + v) = &q,u + 3iq,v and ^n{au) = a ^qu; 

3. if in addition u and v are bounded, then uv G Doni(^Q) with 

dS>n{uv) = V d^nu + u d^nv + 2 {du, dv) dfi. 

Note here that since u, v are in N^''p{Q), it follows that they are well-defined up to sets 
of Cap2-zero; such null sets are not charged by ^uu, ^nv, see Lemma 3.6. 

3.2 Inhomogeneous Dirichlet problem 

In this section we consider the inhomogeneous Dirichlet problem on bounded open domains 
Q such that fi{X \Q) > 0; we assume that a metric space X satisfies a (1, 2)-Poincare 
inequality. More precisely, given two functions / G L^(X, /x) and v G N^''^{X, ^), we wish to 
find u G Dom(Ao) such that 

AfiM = / on f], 

u — V & Nq (f2,/i). 



By definition of An, we interpret ( 1121) in the weak sense, i.e. u is a solution of (1121) if 

u — V E Nq' {Q, /i) and for all (p E Nq' {Q, /i), 

Jn 

As in the classical linear theory, a solution u to (fT2|) can be written as the sum of two 
functions, Uq and Ui, where Uq is harmonic in ^2 such that Uq — v G Nq' (fi,/i) and Ui G 
A^o' (^;/^) is a particular solution to the problem A^ui = f with -Ui G A^q' (^j/^)- 
The function uq is constructed in |KRSj as the minimum of the energy functional 



min / Idul"^ dn- 

eN^''^in,n) Jn 



u-veNg'^{n,ij.) Jn 
For the second part, we use the functional F : N^ (fi,/i) — )■ M given by 

^i'^) ~ o \dufd^+ / uf dfj,, 
2 Jn Jn 

which is the sum of a linear functional and a strictly convex energy. Hence F itself is strictly 
convex. Then, if F has a minimum, it is unique and the minimum is the desired solution 
ui. To prove the existence, it is enough to use the Sobolev inequality, i.e. if m G Aq' (^) A^) 
there exists a constant Cg > such that 

\\u\\2 < Cs\\du\\2. 
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Given that Q is bounded and /i(X \ f2) > 0, the above Sobolev inequahty holds; we refer 
to |HaKo] and |KiSh] for the details. Then, for any u & Nq' (fi,/i) we have that, using the 
inequality ab < ea^ jl + }? jle with a,b,6 > 0, 

F{u) =hdu\\l + fufdf^y hduWl - WfhWuh > hduWl - c^WfUduh 

>(i-f)ll*lli-lll/lll. 

If we fix e < 1/cs, the preceding inequality gives us that F is bounded from below by 

-2-^e-^Cs\\f\\l Therefore, 

m = inf F{u) 

is finite, and in particular, the infimum is a minimum as seen by taking a minimizing sequence 
and applying Mazur's lemma. The minimizing function ui is a weak solution to the desired 
equation, that is 

/ {dui,dip) dfi = — / ffdfi (13) 

Jn Jn 

for all (f E Nq' (fi,/u). From (TT^ it is immediate to see that Ui is the desired solution; in 
addition, if in (TT^ we take (p = Ui, we have the Caccioppoli type estimate 

IM^iilb < c, 

4 Functions of bounded variation and the perimeter 
measure 

The aim of this section is to study some properties of the perimeter measure of a ball in 
metric space. The properties we have in mind are needed in the characterization of a singular 
function which will be constructed in Section [61 

Following [Mrj . the definition of the total variation of a function u G L^{X, /i) is given by 



\Du\{X) = inf < liminf / gu^ d^ : Uj G LipiQ^(X, /i), Uj — )■ -u in Ll^^{X, fi) > . 



(14) 



A function u is said to have bounded variation, that is, u G BV{X,ij,), if \Du\{X) < cxd. 
Moreover, a Borel set E G X with finite measure is said to have finite perimeter if xe G 
BV{X,^). We denote the perimeter measure of E by P{E,X) = \Dxe\{X). 

To each function of bounded variation we associate a Borel regular measure, its total 
variation measure. This measure is defined on every open set A C X using 014p . that is, 

|Dm|(A) = inf < liminf / Qu^ dfi : Uj G Lipi^^{A,^), Uj — ;■ m in LI^^{A, fi) > . 

We extend this measure to act on any Borel set i? C X by the Caratheodory construction 

\Du\{B) = ini {\Du\{A) : A open and B C A} ; 
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for more details on this construction in the metric measure setting see [Mrl Theorem 3.4]. 

An equivalent definition can be also given by way of the Cheeger differentiable structure 
as follows: 

\DcU\{X) = inf < liminf / \duj\ d^ : Uj E Lipi^^{X , fi) , Uj — )■ m in Ll^^{X,fi) 
I J Jx 

and we shall say that u has bounded total Cheeger variation if \Dcu\{X) < oo; a set with 
Cheeger finite perimeter is a Borel set E with finite measure such that \DcXe\{X) < cxd. 

By the results contained in [C], it follows that these two definitions are equivalent, in the 
sense that that u has bounded total variation if and only if it has bounded total Cheeger 
variation. There exists a constant c > 1 such that 

-\Du\{X) < \Dcu\{X) < c\Du\{X). 

Also using the Cheeger differentiable stucture, we have that \Dcu\ defines a finite Radon 
measure; the argument is similar to the case of \Du\ and so we refer to |Mrj for the proof. 

A sequence of Lipschitz functions {uj)j is said to converge in variation to a function 
u G BV{X, fj,) if Uj conveges to u in Ll^^{X, /x) and 



/ gu.dfi^ \Du\{X). 



The preceding definition of total variation does not specify the optimal sequence, i.e. the 
sequence which converges to u in variation. 

It is proved in [Mrl Theorem 3.8] that the discrete convolution gives an approximation 
that is only comparable in variation with the optimal one. Note that for the optimal sequence, 
the vector valued measures dflj = duj dfi, have uniformly bounded total variation. So, up to 
subsequences, they converge to some vector valued finite measure /2oo. 

Remark 4.1. A sequence {uj)j converging to u in variation is optimal not only for the 
variation in X, but also for the variation in all open subsets A with \Du\{dA) = 0. In fact, 
by definition, we have that 

|Z^m|(A) < liminf / g^ dfx, 

J^°° J A ' 

but also that 



\Du\{X\A) < limsup / _5f„^ dfi < limsup I / gu^ d^ - I gu^ djj, 

i-s>oo Jx\A j^oo \Jx J A 

= |L)m|(X)- liminf [ gu dfi < \Du\{X\A). 
i~^°° J A ' 

The preceding inequalities are indeed equalities if \Du\{dA) = and so \Du\{X \ A) = 
\Du\{X \ A). Hence the following two limits exist 

lim / gu,dfi= \Du\{X\A), 
•?'^°° Jx\A 

and 

lim / gu^dfi= \Du\{A). (15) 

^^°°Ja 
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An important tool in the theory of functions of bounded variation is the coarea formula. 
The version we work with in the present paper is a direct consequence of |Mrl Proposition 4.2]. 
For any u G BV{X, fi) and any Borel measurable function / : X — t- M, the following identities 
hold 

/ fd\Du\ =11 fix)d\DxEMx)dt (16) 

Jx JR Jx 



and 



fd\Dcu\= / / f{x)d\D^XEt\ix)dt, 
X Jr Jx 



where Et = {u > t} , t E R, is the super- level set of u. We point out that in these formulae, 
due to the fact that the measures \Du\ and \Dcu\ are not absolutely continuous with respect to 
the measure /i, it is important to consider the function / and not an equivalent representative. 
Since the perimeter measure does not charge sets with zero 1-capacity, we can modify the 
function / on such negligible sets. If u is Lipschitz, f ll6p can be written as follows 



fgudn= / / f d\DxEt\{x)dt. 

jR Jx 



This follows by an argument contained in |Cat Theorem 6.2.2] and summarized in the fol- 
lowing proposition. We will provide a proof here for the reader's convenience. 

Proposition 4.2. Let u G Lip(X). Then the total variation measure d\Du\ is given by Qudfi. 

Proof. Lipschitz continuity of u implies that \Du\ is absolutely continuous with respect to 
fi with density given by some function Gu- To see this, note that Uj = u is a. possible 
competitor in the definition of \Du\, and so for every open set A G X (and hence for every 
set y4 C X) we have 

\Du\{A) < / gudfi. 

J A 

Therefore the density function Gu < gu /i-a.e. 

To prove the equality it suffices to prove that the function G„ is an upper gradient of u. 
We take a sequence {uj)j of Lipschitz functions converging to u in V-{X,jj) and with 



lim / gu.dji = \Du\{X). 



3^°° Jx 

The sequence of measures gu dfi is bounded, and so, up to a subsequence, it converges weakly 
to a measure /loo that is still absolutely continuous with respect to fi, i.e. dfi^o = goo dfi. 
To see that fi^o is absolutely continuous with respect to fi, it suffices to show that whenever 
E C X is compact with fi{E) = 0, we have fioo{E) = 0. To this end, we note that because E 
is compact, for every e > we can cover E with a finite number of balls B^ with ^{dBD = 
so that the open set A^ = Uj-Bf contains E and is such that ^{Air) < e and jjii^dA^) = 0. 
It follows that iDuKdAgr) = because of the absolute continuity of \Du\ established above. 
Therefore, by Remark 14. II we have 



lim / 5f„, djj. = \Du\{A^) < gudjj. < Le, 
^■^~ J A, J A, 



which imphes that fiooi^e) < Le. It follows that /ioo(-E) = 0. For more general sets E G X 
with /i(-E) = 0, there is a Borel set £"0, containing E, such that /i(-Eo) = 0. Because /ioo is a 
Borel measure, /ioo(-Eo) is the supremum of all /ioo(-ft'), the supremum taken over all compact 
sets K C Eq. Given that fi{K) = and so fiooiK) = 0, it follows that /ioo(-Eo) = and so 
lioo{,E) = 0. Then gu^ converges to goo weakly in L^(X, /i). 

To summarize, we have Uj — )■ m in L^{X) and gu. — > goo weakly in L^{X). Now, an 
invocation of the Mazur lemma, together with [KaSj Lemma 3.1] shows that goo is a weak 
upper gradient of u; it follows that g^ < goo a-e. in X. 

Given that balls have finite /i-measure, when x E X, for almost every r > we have 
fi{dBr{x)) = 0. For such r > 0, by ([I5]) with A = 5,(x), 

\Du\{Br{x)) < / good^ <\immi / (7^^. (i/x = |DM|(i?r(a;)) = / Gudfj,, 

and so we get goo{x) = Gu{x) for x G X that are Lebesuge points for both Gu and goo- It 
follows that Gu < gu < goo = Gu a.e. in X, from which the claim follows. 

D 

From now on, we assume that X is also a geodesic space. This is not an overly restrictive 
assumption, since X, by the virtue of supporting a Poincare inequality and being complete, 
is a quasiconvex space. It follows that in our setting, a bi-Lipschitz change in the metric 
does result in a geodesic space. 

Let us fix a point Xq G X. Since X is assumed to be a geodesic space, by the results 
in [C], the function Uxo{x) = d{x,Xo) is Lipschitz with Lip^u^o) = gu^ = 1- Moreover, we 
may write -Bt(xo) = {m^o < t} and so by the coarea formula fITBl) we obtain for any positive 
r > that 

/ PiBtixo),X)dt = ii{Brixo))<oo. 
Jo 

Thus the map t 1— ;■ P{Bt{xQ),X) is a measurable locally integrable function. This implies 

that for almost every r > 0, 

P{Br{xo),X) = \im- P{Bt{xQ),X)dt = \im- P{Bt{xo),X) dt 



' r—e 
"r+e 



= lim;i / P(5i(xo),X)dt. 
e^o 2e J^_, 

In particular, for almost every r > the perimeter measure coincides with the Minkowski 
coiiijGiiij 

For a ball Bri^Xo) satisfying (TT7|) . we can consider the sequence of functions (M£)e>o, where 



r — d{xo,x) 1 „ 1 . f 1 



Me(a;) = max ■{ min <( ^ — , 1 > , > = min < -d{x, X \ Br{xo)), 1 f . (18) 



For a such function u^, we have that (^^^ = -XBr{xo)\Br-e{xo) and 

gu,djj,= - djj,-^ P{Br{xo),X), 

JX ^ J Br(xo)\Br^e(xo) 



' Br{xo)\Br~e{xo) 
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that is, the sequence u^ converges to XBr{xo) in variation. This also means that the sequence 
of vector valued measures {Idueld^)^ is equibounded 



\DcUi;\{X) = / \dui;\ d^ < c < oo 
Jx 



for some positive constant c. Therefore there exists a subsequence Sj — )■ such that, setting 
Uj = Me . , the sequence of vector- valued measures duj dfi is weakly convergent to some vector 
valued measure /Too- This measure is absolutely continuous with respect to both \DcXBrixo)\ 
and \DxBr{xo)\- Indeed, if, for instance, \DcXBrixo)\iE) = 0, where E C X is compact, 
we can find for every e > an open set A^ D E such that \DcXBr{xo)\{dAs) = and 
\DcXBr{xo)\{As) < e. Reasoning as in the proof of Proposition 14.21 we may conclude that 

|/ioo|(^) < |/ioo|(^e) < liminf / \duj\ dfi = \DcXb,.{xo)\{A) < £■ 
We may hence write 

Atoo = ^xo,r\^cXBr{xo)\ = 0'xo,r\DXBr{xo)\ 

for some vector-valued \DcXBr{xo)\-'^^^^'^T^^^^^ function z/^^^r and |DxBr(a;o)l"™^^^^^^t)le func- 
tion (Txo.r- In particular, the function ax^.r plays the role of the normal vector at the boundary 
of Br{xo). In this context, it is not clear from the definition if it is a unit vector. For the 
sake of simplicity, if no confusion may arise, we simply denote the functions i'xo,r and a^Q^r 
by z/ and a, respectively. 

We can summarize the previous construction in the following definition. 

Definition 4.3. We shall call a ball Br{xo) regular if the equation (TT71) is valid and if there 
exists a sequence Ej — )■ such that for the sequence of functions Uj = u^^ , referred to as an 
optimal sequence and defined in flTSl) . the following hold true: 

(1) gu^ dfi converges weakly to d\DxB,.{xo)\', 

(2) duj dfi converges weakly to a d\DxBrixo)\ fo^ some |DxBr(a;o) [-measurable vector- valued 
function a = <7xo,r- 

Almost every ball is regular in the sense that for every xq G X and for almost every r > 
the ball Br{xo) is regular. Howeover, the vector a is not a priori unique and it is not clear 
whether it depends on the sequence Ej we consider. 

The given notion of regularity relates to interior regularity of a ball. One can also consider 
the notions of outer and two-sided regularity and obtain that for almost every radius r > 
the ball Br{xo) has inner, outer, and two-sided regularity. 

5 Divergence measures and generalized Gauss— Green 
formulas 

Here we consider divergence-measure fields, i.e. a class of vector fields F : X —> M.^ belonging 
to the space L'^{Q,'R^, fi) and for which divF is a measure. In the metric space framework 
of the present paper, we generalize some results obtained by Ziemer in [Z]. 
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Previously, Thompson and Thompson in |TT j constructed a divergence form in the setting 
of Minkowski spaces, and they proved a Minkowski space analogue of the Gauss-Green 
theorem. 

The aim of this section is to study the operator div on Lip^(r2) with values in the space 
of measures, that is, we want to define for F G L^(r2,]R'^,yu) the distribution 

(m, div F) :=- I {F, du) dfi, (19) 

Jn 

for u G Lip^(r2). In the following we adopt the notation from |CTZj and [Z]. 

Definition 5.1. We say that F G L^(r2,M'',yu) is in the class PA^^(i7) if there is a signed 
finite Radon measure, denoted by divF G ^^(fi), on fi such that 

/ uddivF = - {F,du)dfi (20) 

Jq Jn 

for all u G Lip^(r2). 



Remark 5.2. By an argument similar to the proof of Lemma 13.61 we can prove that if 
F G VAi'^^Vl), then the measure divF does not charge sets with zero 2-capacity. Therefore, 
condition f l20|) can be extended to any u G Nq' (f2,/i) n L°°{fl,fi). 

Note that when F G L'^{n,R'', fi) the operator Tp : A/'o'^(fi,/i) -^ M given by 

Tp{u) = / (F, du) d^ 
Jn 

is a bounded linear operator on the Hilbert space (A'"q' {Q,fi),Si). Therefore, by the Riesz 
representation theorem, there exists a function f G A^g' (^)/^) such that whenever u G 
No''^{n,fj,), Tpiu) = £iiv,u). Hence, if F G VM^{n), then for all u G A/'o'^(f],/i) we obtain 



that is, 



{F, du) dfi = / uv dfi+ / {du, dv) dfi, 
n Jn Jn 



uddivF+ uvdfi = — {du,dv)dfi. 
n Jn Jn 



It follows that V G Dom(^Q) with 

d^QV = —V dji — d div F. 
This proves the following lemma. 

Lemma 5.3. Given a domain Vt <Z X , a map F G L^(r2,M'^,/i) is in the class VM^iyt) if 
and only if there exists v G Dom(^Q) such that 

dS>QV = —V djji — d div F 
in the sense of distributions on Nq' {Q,fi). 
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We can also state the following simple properties of the divergence measure. 

Lemma 5.4. Let F G VM'^{^). Then spt(divF) C spt(F). Moreover, if v e Dom(^Q), 
then dv G VM'^{^) with div dv = ^qV. 

Proof. The first statement follows by considering A = f2 \ spt(-F), so we have that 



div F\{A) = sup <^ / ipddivF : if e Lip^(y4), ||v9||oo < 1 
= sup < / (F, dip) dfi: (p e Lip^(y4), \\(p\\oo < 1 



0. 



For the second part, let v G Dom(^n)- Then there exists a signed finite Radon measure 
^nv G ^b{X) such that 

{dv, du) d^= — I u d^QV 



for all u G Nq' (fi,yu) fl L°°(i7,/i). From this and by Remark 15.21 we may conclude that 
dv G VAi'^{Vl) with divfii; = ^nv, and the claim follows. D 

We now state the following two propositions on the Gauss-Green type integration by parts 
formula for vector fields in VM°°{n), that is for vector fields F in L~(l], M^ /i) n VM^{n). 

Proposition 5.5. Let F G VM°°{n) n C{n,m'') and Br{xo) d Vt he a regular hall. The 
following Gauss-Green formula 



fddwF+ / {F,df)dfi=- / f{F,a^^^r) d\DxBr{xo)\, 

Brixo) J Br(xo) Jn 

holds for all f G Lip^(r2). If the support of F is disjoint from dBr{xo), then the requirement 
that F is continuous can he removed. 

Proof. We can consider an optimal sequence of locally Lipschitz functions {uj)j converging 
to XBr{xo) ill variation as in flTSl) . Then we have by the Leibniz rule that 

/ UjfddwF = — / {F,d{ujf)) dfi = — / Uj{F,df) dfi— / f{F,duj)dfi. 
Jn Jn Jn Jn 



We notice that 



UjfddivF- / fddivF 



< 



divF 



{Br{x^)\Br-e,{xo)), (21) 



and that the right-hand side of fl2Tl) tends to as j — t- oo. Thus we may conclude that 



lim / UjfddivF= / fddivF. 
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Also, by the fact that both F and df are in L°°(fi), by an apphcation of the Lebesgue 
dominated convergence theorem we obtain 



hm / Uj (F, df) d^= I (F, df) d^. 

^^°°Jn JBr-ixn) 



We also have, due to the continuity of F, that 

lim / f{F,duj)dfi= / f{F,a^g^r)d\DxBr{xo)\, 
and so the proof is completed. D 



Remark 5.6. We point out that property (1211) is a consequence of the choice of an optimal 
sequence {uj)j to be an inner approximation of the characteristic function XBr{xo)- If we 
chose, for instance, an outer approximation, then the preceding integration by parts formula 
would be as follows 

/ /cidivF+ [ {F,df)dii=- [ f{F,a^^,r)d\DxBrixo)l 

JBrixo) JBrixo) Jn 

for all / G Lipg(f2), where axo,r is the density of the vector- valued measure obtained as a 
weak limit by way of the gradients of this new sequence as in Definition I4.3[ 

We also point out that the previous proposition can be extend to more general sets 
E G Q with finite perimeter whenever a Minkowski content characterization of the perimeter, 
analogous to flT7|) . holds. In this case, the boundary of E has to be considered as the essential, 
or the measure-theoretic, boundary of E, i.e. the set of all points at which the density of E 
is neither nor 1. 

We prove the following main theorem of this section, which is a generalization of Propo- 
sition 15.51 without requiring continuity of the vector field. This theorem should be thought 
of as the generalization of the Gauss-Green theorem of the Euclidean setting. 

Theorem 5.7. Let F e VM'^{^) and let Br{xo) <zVL he a regular hall. Then the following 
extended Gauss-Green formula 

I fddivF+ [ {F,df)dfx= [ f{F-u)- d\DxBrixo)l (22) 

JBrixo) JBrixo) Jn 

holds for all f G N^'^{Q,fi) fl L°°{Q,fi), where {F ■ v)~^Brix ) '^^ ^^^ interior normal trace of 
F on dBr{xo). 

Proof. We use the optimal sequence {uj)j defined in (TT8|) . Then, as in the proof of Proposi- 
tion 15. 5[ by the definition of divF (Definition 15. ip and the Lebesgue dominated convergence 
theorem 

lim / f{F,duj)dfi= lim I / {F,d{ujf)) dfi — / Uj{F,df) dfi 
:>^°° Jn J'^°° \Jn Jn 



— lim I / UjfddivF+ / Uj{F,df) dfi] 
J'^°° \Jn Jn J 



fddivF- / {F,df)dn. 

Br ixo) JBrixo) 
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For the sequence {Lj)j of operators given by 

Ljif) ■= / f{F,duj)d^, 
Jn 

we have that \Lj{f)\ < C||F||oo||/||oo, where the positive constant C is given by 

C = sup / \duj\ dfi < oo. 

Indeed, C is finite since the ball Br{xo) has finite perimeter. In particular, C is independent 
of both / and F and so, by the above argument, the operator 

L{f) := hm L,(/) 

is bounded over Lip^(i7) and admits an extension to Cc{^). This in turn implies that there 
exists a measure z/ G ^^(fi) such that for any / G Cc{^) 



L{f)= [ fdu. 
Jn 



The measure z/ is concentrated on dBr{xo)] in fact take any compact set K such that 
K n dBr{xo) = 0, an open set A D K such that dist {A, d Br {xq)) > 0, and take Ej < 
dist(A, dBr{xo)). Then, since spt{duj) fl A = 0, we obtain for any / G Lip^(y4), 



f du = lim / f{F, duj) dfx = 0, 
Jn J-^°° Jn 



that is Iz^l(-ft') = l'^l(^) = 0. This property extends to any Borel set E such that E fl 
dBr{xo) = since 

\iy\{E) = sup \iy\{K) = 0. 



KCE 



Also z/ can be seen to be absolutely continuous with respect to \DcXBrixo)\] indeed, if E 
is a Borel set such that \DcXBrixo)\{E) = 0, then there exists an open set A^ such that 
\DcXBr{xo)\{^e) < ^- Fix a compact set i^ C -E and an open set A D K such that A C A^. 
Then, for any / G Lip^(A) with ||/||oo < 1 we have that 

/ f{F,duj)dfj, < limsup||F||oo / \duj\ dn < \\F\\^\DcXBr(xo)\{^) < ^Wf^Woc, 

that is |z^|(A) < e. Therefore, since e is arbitrary, jz^K-ft') = 0. Finally, by taking the 
supremum over K G E, we obtain that Iz/K-E) = 0, and hence z/ is absolutely continuous 
with respect to \DcXBr{xo)\- 

To conclude, there exists (F ■ ^)QBr(xo) ^ ^^i\^XBr{xo)\) such that 



Lif) = - [ fiF-^)-aB.ixo)d\DXB.ix,)\ 



— * 

This map defines, in the metric setting, the interior normal trace of F on dBr{xo), and the 
integration by parts formula (l22l) holds. D 
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Remark 5.8. The term interior normal trace can be justified by the following facts. If 
F G VM^^iVt) n C((],M^), then by Proposition [53] we get that 

In addition, also when F is not continuous, recalling that with Uxq{x) = d{x,XQ), 

duj{x) = -—du^^{x)xBr{xo)\Br-,Axo){^)^ 

and by using the coarea formula (IT6|) . we can write 

1 



/ (F, duj)dij, = / / (F, du^o) dfi 

= / / f {F,du^^) d\DxBtixo)\dt. 

Therefore, we have obtained that 

/ /(^■'^)ai3,(xo)^I^^B'-(^o)l = - lim — / / f{F,du^^)d\DxBt{xo)\dt, 
Jn ^ ' J^°° £j Jr-sj Jn 

which gives meaning to the following equality in terms of the trace 

/ fi^-'^)dBAxo)^\^^Brixo)\ = - / f{F,du^,,)d\DxBr{xo)\, 

and to the fact that the vector du^o defines in a weak sense the normal vector axo^r to dBr{xo). 

Remark 5.9. Observe that in the proof of Proposition 15. 71 we have used a particular optimal 
sequence. It turns out, nevertheless, that the interior normal trace (F ■ J-')^Br{x ) "^^^^ '^^^ 
depend on this particular choice. This fact is a direct consequence of equation ( !22l) . since 
then formula 

[ fiF-^rdBAxo)d\DXBAxo)\= I fddivF+ f {F,df)dfx 

uniquely identifies the values of (F ■ J-')^^^/^ y 

Remark 5.10. By ^ Theorem 5.3] (see also |AMP] ) . formula (122|) can also be written by 

/ fddivF+ [ {F,df)dfi= [ f{F-v)- d.,,rdS\ 

JBrixo) JBAxo) Jd*Br{xo) 

where d*Br{xo) is the essential boundary of Br{xo), S^ is the spherical Hausdorff measure 
defined using the Caratheodory construction based on the gauge function 



hiB,) 



KB,) 



and "(^xo^r : X — )■ [c, Cd] is a Borel function depending, in general, on the ball Br{xo), and c is 
a positive constant and q the doubling constant of /i. 
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6 Harmonicity and the mean value property 



In this section, we shall follow the approach of |HS] and construct, for any regular ball 
Br{xo) C X and any x G Br{xo) the Green function on Br{xo) with singularity at x, that is 
an extended real-valued function G{x) = G^ /^^Jx) such that 

1. G is strictly positive and harmonic in Br{xo) \ {x}; 

2. Ge N''\X \ B,{x)) for any e > and G\^^s^^,^^) = 0; 

3. for every y G dB^^xo) 



4. G is singular at x; that is 



limG(x) = 0; 



liniG(x) = cxd; 



5. for all < a < 6, 

Cap2({2; G BriyXo) : G{x) > b}, {x G Br{xo) : G{x) > a}) 



b — a 



In |HS] the authors constructed the Green function of a relatively compact domain with 
the aforementioned properties in metric measure spaces; we refer also to [Hj and |DGM] . We 
can state the existence and main properties of the Green function in the following theorem. 
We assume that X supports a (1, 2)-Poincare inequality. 

Theorem 6.1. Let Q G X be a relatively compact domain. Then there exists the Green 
function G = Gf^ with singularity at x G Q. In addition, dG G L^(X \ ^^(x)) for any e > 
and 

where Uq is a positive Radon measure in the dual Nq' {X \ B^^x))* concentrated on dVt. 
Moreover, G admits the measure-valued Laplace operator 

S)xG = 4 - vn. 
in the sense that for any v G A^^'^(X) continuous at x, then 

{dG,dv)dfi= / vdv^—v{x). 
X Jdn 

Proof. We refer to |HS] for the details on the construction of G. We sketch the main steps 
needed in the definition. We find a harmonic function on Q\Bi; (x) 

G -' 

' Cap2(5,^.(x),fi)' 

where B^.{x) is a regular ball, Sj \ 0, Sj < dist(x,(9fi), and Vj is the potential of B^.{x) 
with respect to Q; that is Vj G A^^'^(X) is harmonic in Q\Bs{x), Vj = OonX\Q and Vj = 1 
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on Be {x). It is then shown that, up to subsequences, the functions {Gj)j converge locally 
uniformly in X \ {x} to a function G. The limit function G has the desired properties of a 
Green function. 

Let us fix a positive sequence (Mj)j>o such that Mi /^ oo, and the truncations 

TiG:=min{G',Mi}. 

There exists a sequence r^ \ of radii such that 

Ei cBrXx), 

where we have written E^ = {x & Q : G{x) > Mi}; and we may consider the case in which 
Tj < e. Then TjG = G on X \ B^i^x) and TjG is subharmonic in X \ B^ix). By |BMSj 
(we refer also to [M] for a detailed description in the Euclidean case) there exists a positive 
Radon measure u^ in the dual Xq' (X \ Be{x))* such that for all v G Lip^(X \ -Be(x)) we 
have 



/ {dG, dv) dfi = V duQ . 

Jx\Be{x) Jx\B^(x) 

If f G Lip^(X \ n), the fact that G = on X \n implies dG = on X \ fi, and then 

/ V dv^ = / vdi'Q= / {dG,dv) dfi = 0. 

Jx\B^{x) Jx\n Jx\n 

On the other hand, the harmonicity of G in fi \ -Be(x) implies that ii v E Lip^(fi \ Be{x)), 

V du^ = V du^ = / {dG, dv) dfi = 0. 

x) Jn\B^{x) Jn\B^{x) 



then 

IX\B,(a 

Hence the measure z/^ is concentrated on dfi. 

Analogously, since TiG is superharmonic in Q there exists a positive Radon measure 
i/f G Xo'^(fi)* such that for all v G Lip^(fi) 



/ {dTiG, dv) dii=- I 



V uu 

Q\Br; (x) 



vdvf. 



The measures Vi are supported in Br^{x); indeed, since TiG = G on f2\i?j..(x) it is harmonic. 
Hence, ii v E Lip^(fi \ Br..{x)), 

/ v dvf = V dvf = I {dG, dv) dfi = 0. 

Jn J^\Br^ (x) Jn\Bri (x) 

Following the argument of Serrin [SeJ Lemma 1 and Theorem 3], there exists A G M such 
that if f G Lip^(i7) is equal to 1 in a neighborhood of x, then 

/ {dG, dv) dfi = A. 
Jn 
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Indeed, if fi,f2 G Lip^(r2) are two functions that are equal to 1 in a neighborhood of x, the 
difference f = Wi— ^2 belongs to Lip^(r2\{x}); hence, the harmonicity of G in 0\{a;} implies 
that 

/ {dG, dvi) dfi — {dG, rff 2) dfi = {dG, dv) djj, = 0. 
Jn Jn Jn 

In particular, if t> G Lip^(r2) is a function such that f = 1 on Brj^(x), then 

ly^CBrXx)) = V duf = - {dTiG, dv) dji = - I {dG, dv) dji = -A. 
Jn Jn Jn 

This argument implies that A G M is negative and the measures uf^ are equibounded in 
^5(^2) . Thus, up to subsequences, uf converges weakly to A^^. 

To summarize, we have proved that the sequence of the measure- valued Laplace operators 

^T,G = uf - uE 

admits a convergent subsequence ^Ti^G, defining the measure- valued Laplace operator 



:G = hm ^T,^G = \6, - z/g. 



The fact that the limit measure is uniquely determined implies that for any sequence Mj /^ 
00, the measures S^TiG converge and the limit measure is \5x — i'q- 

Let us show that A = —1. Let us consider the set i? = {x G f2 : G{x) > 1} and a 
function v G Lip^(r2) such that f = 1 on ii^. Since x is an interior point of E, we have 



A = Xv{x) = — {dv, dG) dji. 

Jn\E 



On the other hand, the map f = {G — v)xn\E belongs to Nq' (fi \ E) and then 

= / {df, dG) dfi = |(iGp d^ — {dv, dG) d^. 

Jn Jn\E Jn\E 

These properties of G imply that G is the potential of E with respect to Vt, that is 



Jn\E 



in\E 

We may hence conclude that A = — 1. 
Finally, we point out that the identity 

{dv, dG) dfi = V du^ — v{x) 
Jan 

is valid for functions v G N^'^{X) that are constant in a neighborhood of x, but it can be 
generalized to functions v G N^'^{X) that are continuous at x. This is a simple consequence 
of the limit 

an *^°°ian Jsrix) 



■ lim / {dv, dTiG) dfx = {dv, dG) djj,. 
*^°° Jx Jx 



D 
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Remark 6.2. Let us consider the (first) Heisenberg group EI with the geodesic distance. In 
this case, the natural differential structure is given by the horizontal bundle and the Laplace 
operator is just the horizontal Laplace operator. In this setting, we can use all the results 
of the preceding section and obtain the representation of the measure z/*^ in terms of the 
perimeter measure. Notice that a ball Br{xo) in HI satisfies a ball condition as in |AKSZ[ 
Definition 2.1] at its boundary except at two points; a finite collection of points is negligible. 
Hence, if G = G^^,^ -, is the Green function on Br{xo) with singularity at x, then whenever 
X is a boundary point of Br{xo) satisfying the ball condition, 

^(G, X, g) := sup G - sup G < Cg, (23) 

where < p < d{x,x)/2 and G is a positive constant that does not depend on x, x, or p. 
It follows from a covering argument together with (l23l) and |BMSl Lemma 4.8] that z/^ /^ -, 
is absolutely continuous with respect to the perimeter measure \DxBr(xo)\- Moreover, there 
exists a function {}g G L^{X, \Dxb\) such that du^ = ■dcdlDxBl- The function ^g comes 
from the Radon-Nikodym theorem. 

We give a characterization of harmonic functions via a mean value type property with 
respect to boundary measures. 

Theorem 6.3. Let u G N^''^{VL,ij), then the following hold: 

(1) Let u he harmonic in Vt. Then for every regular hall B^.{xq) C VL and x G Bj.{xq) 

u{x)= I wrf^f.(.o); (24) 

JdBr(xo) 

(2) If for every regular hall Br{xQ) C Q and any x G Br{xQ), u satisfies the mean value 
property ( l24l) . then u is harmonic in Vt. 



An analogous characterization holds true for suh- and superharmonic functions. Let u G 
N^''^{Q,fi) then the following are equivalent: 

(3) Let u he suhharmonic (superharmonic) in Vt. Then for every regular hall Br{xo) C Q 
and X E B 

JdBrixo) V JaBr(xo) 

(4) If for any regular hall Br{xo) and any x G B^{xq) 

u{x)< I Mc?z/f( ), \u{x)> I udu^ 



dBr{xo) \ JdBr(xo) 



then u is suhharmonic (superharmonic) . 
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Proof. Suppose that u is harmonic. Then u G A^^'^(r2,/i) fl L^^{Q) and we can apply Theo- 
rem [6111 We obtain for any regular ball Br{xo) and x e Sr(^o) 

= / {du,dG%) dfi= - ud&xG% = -u{x) + / wrfz/f , 

Jx JX JdBr(xo) 



which gives the condition (1). 

On the other hand, if u is continuous, if we fix a regular ball B = Br{xQ), we can consider 
the harmonic function Hu generated by u on B, that is the solution of the problem 



min < / \dv\'^ d^ w — u E Nq' {B, fi) > . 



Hence Hu is harmonic in B and satisfies the mean value property, that is for any x E B 

Hu{x) = f HJu^. (25) 

JdBrixo) 

The conclusion follows from continuity of u since 

lim HJx) = u(y), 

B3x^yedB ^ ^ ^^^' 

and then by ( l25l) . Hu = u on B. For a general u G N^''^{Q, fi), we can find a continuous 
function u^ such that u = u^ on a. set of capacity less than e and such that \\u — M£||i,2 < ^\ 
then via a standard approximation argument, we can conclude. 

The same line of reasoning carries out in the case of sub- and superharmonic functions. 

D 

Remark 6.4. It was proved in |BBSj that the harmonic extension of a function u G 
A^^'^(r2,yu) on a ball B (Z Vt can be expressed in terms of harmonic measures Vx with sin- 
gularity at a; G -B; by this we mean that if (/J G C{dB), then in jBBS| Theorem 5.1], its 
harmonic extension is given by 

H^{x) = / (pdus;. 

JdB 

If we move x G Br{xo), it is possible to see that the measures u^ are mutually equivalent; 
in particular, if we take xq and x G Br{xo) \ {xq}, we have that Ux is absolutely continuous 
with respect to Ux^ and its density P{x, ■) is called the Poisson kernel. In other terms, the 
Poisson kernel is defined as 

P{x,x) = -— ^(x). 
dux„ 

In |BBSj . Ux was not explicitly identified. Nevertheless, from the results contained in the 
previous sections, we are able to identify this measure as the outward normal derivative z/^ 
of the Green function. 



Example 6.5. In Example 13.51 if we take fi = -Bi(O), the unit ball, then all balls except 
-Bi(O) are regular. This is due to the fact that the perimeter of -Bi(O) has weight 1, that is 
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I-DX-bI = ^^ ^\-dB. However, if we consider the optimal sequence {uj)j defined in (fT8|) we 
have that 

Nevertheless, the measure z/^ /q-, can still be characterized as a perimeter measure, but with 

dz/f^(o) = 2{VG ■ z/Bi(o))rfH""'L9fii(0) = 2{VG ■ z/B,(o))t/|/^XB,(o)|. 

On the other hand, if we take any other ball B C -Bi(O), it is regular and in this case 
W-\dB n OBiiQ)) = 0. Note also that if K^-^B n dBi{0)) > 0, then since the Green 
function is harmonic in B except the singular point x, we have that VG ■ i^Bi{o) = and then 

du'i = (1 + Xb,(o))(VG ■ UB)d'H''-'LdB = (VG • UB)d\DxB\ 

On the other hand, if we take fi = M" \ -Bi(O), then every ball is regular. This is due to the 
fact that in this paper regularity is a notion of inner regularity. If one changes the notion to 
outer regularity or to two-sided regularity, then things change. 
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